Quantum Effects in the Mechanical Properties of Suspended Nanomechanical Systems 
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We explore the quantum aspects of an elastic bar supported at both ends and subject to 
compression. If strain rather than stress is held fixed, the system remains stable beyond the 
buckling instability, supporting two potential minima. The classical equilibrium transverse 
displacement is analogous to a Ginsburg-Landau order parameter, with strain playing the role 
of temperature. We calculate the quantum fluctuations about the classical value as a function of 
strain. Excitation energies and quantum fluctuation amplitudes are compared for silicon beams 
and carbon nanotubes. 
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The continuing drive towards semiconductor device 
miniaturization and integration has resulted in fabrica- 
tion and micromachining technologies that are capable of 
producing artificial structures with features approaching 
the ten nanometer length scale. To go beyond this scale, 
naturally occurring and chemically organized structures 
are receiving much attention. The availability of these 
top-down and bottom-up nanofabrication capabilities has 
initiated the new area of nanomechanics ||, |[ ^, ||] 
in which ultra small mechanical systems are used to ex- 
plore both fundamental and applied phenomena. Re- 
cently, two reports have appeared on two-state nanome- 
chanical systems. In one crossed carbon nanotubes 
were suspended between supports and the suspended el- 
ement was electrostatically flexed between two states. In 
the second , it was proposed to use an electrostatically 
flexed cantilever to explore the possibility of tunneling in 
a nanomechanical system. 

In this Letter we discuss quantum effects in a two-state 
mechanical system that has a tunable, symmetric poten- 
tial function. This mechanical system has analogies to 
the superconducting interference device in which the first 
observation of a coherent superposition of macroscopi- 
cally distinct states was recently reported l^. Specifi- 
cally, we consider a suspended elastic bar under longitu- 
dinal compression. The compression is used to adjust the 
potential energy for transverse displacements from the 
harmonic to the double-well regime, as shown in Fig. |^, 
with strain playing a role analogous to temperature in a 
Ginzburg-Landau system. As the compressional strain is 
increased to the buckling instability ||] , the frequency of 
the fundamental vibrational mode drops continuously to 
zero. By controlling the separation between the ends of 
the bar, i.e. fixing the strain, the system remains stable 
beyond the instability and develops a double well po- 
tential for the transverse motion. Since both the well 
depth and asymmetry are tunable, a variety of quantum 



phenomena may be explored, including zero-point fiuctu- 
ations, tunneling, and coherent superpositions of macro- 
scopically distinct states. In the latter two cases, the sys- 
tem may provide a mechanical realization (at least in the- 
ory) of models studied in Refs. Q and respectively. 
We have applied the model to suspended silicon beams 
and carbon nanotubes, and show that in both cases the 
quantum fiuctuations in position approach 0.1 A; an or- 
der of magnitude greater than the relaxed values. We 
argue that compressing a carbon nanotube at low tem- 
perature may cause a crossover between the quantum and 
thermal fluctuation regimes. Further, we suggest that 
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FIG. 1: Potential energy V as a function of the fundamental 
mode displacement Y. The shape of the potential energy 
is harmonic for e > Ec, quartic for e = Ec = critical strain 
{ec < 0), and a double- well for e < Ec. 
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tunneling in this mechanical system may be difficult to 
achieve because it will require precise regulation of the 
applied strain. 

As a starting point, we consider the normal modes 
and associated quantum properties of an elastic rectan- 
gular bar of length I, width w and thickness d satisfying 
{I >> w > d), supported at both ends without strain. 
We suppose that d is smaller than w so that transverse 
displacements only occur in the "d" direction. The 

equation of motion for small displacements is jl^ 



(1) 



where fi — m/l is the mass per unit length and J- is the 
linear modulus (energy per unit length) of the bar. T is 
related to the elastic modulus Q of the material by = 
Qwd. The bending moment n is given by = for 
a bar of rectangular cross section. 

The normal modes of the bar, y{x,t) = 
Unix) exp{±iujnt) , have the general spatial depen- 
dences 



yn{x) = a„ cos knX + bn cosh fc„x, 
Unix) = an sin knX + bn sinh fc„a;. 



71=1,3,5,. ..(2) 
n = 2,4,6,...(3) 



with an/bn and wavenumbers fc„ fixed by the boundary 
conditions. The boundary conditions for hinged end- 
points, y{±l/2) = = y"{±l/2), lead to normal modes 
with bn — 0, kn — mr/l, and (angular) frequencies 



UJ„ 




(hinged b.c). (4) 



Clamped endpoints have boundary conditions y{±l/2) — 
= ?/'(±Z/2); both a„ and 6„ are nonzero, and the normal 
mode frequencies are given to good approximation by 
replacing n with (n -|- |) in Eq. (^. Because of their 
simplicity, we shall refer to hinged boundary conditions 
when specific results are quoted. 

The mean square displacement of the bar (both zero- 
point and thermal) is an incoherent superposition of con- 
tributions from each normal mode. At the center, only 
even-parity modes contribute, so that 



odd 



2m*iUn 



[l + 2fihL,n/kT)], (5) 



where f{x) = l/(e^ — 1) is the thermal excitation number 
of the n-th mode. For hinged boundary conditions, m* = 
m/2 exactly, whereas in the clamped case m* is slightly 
smaller and weakly mode dependent. In either case the 
fundamental mode is responsible for more than half the 
total mean square displacement. 

Longitudinal compression of the bar moves the fre- 
quencies downward, with a corresponding increase in the 
zero-point motion. Compressive or tensile strain con- 
tributes the "elastic" potential energy Ve = {J- /2lo){lt — 



lof, where k = J dxyjl + {y'Y + dx{y'Y is the 
total (dynamic) length of the bar, I is the endpoint sep- 
aration, and lo is the unstressed equilibrium length. We 
subtract the static contribution {!F/2lo){l — loY since it 
contributes nothing to the dynamics, and add the "bend- 
ing" contribution Vb ^ J dx{y"y to get 



V[y{x)] dx{TK'{yy+Ts{yr) + ^{ f dxiy'f) 



T 

_ (6) 

where e = {I — lo)/lo is the strain, positive if ten- 
sile and negative if compressive. From the Lagrangian, 
L[y{x,t)] — {fi/2) J dx (y)^ — V[y{x)], we find the equa- 
tion of motion. 



(4) 



^ey"-lT 



dx'[y'{x')r]y" = 0, (7) 



which generalizes Eq. (^ . The third term represents the 
tension induced by externally-imposed stretching, and 
the anharmonic fourth term is the enhancement of ten- 
sion due to dynamic stretching; this term arises from the 
geometry of the system. 

In the harmonic regime where where the fourth term 
can be neglected, the normal mode frequencies under 
hinged boundary conditions are given by 



(i+s(—y), (8) 



where ujn are the relaxed {I = lo) frequencies of Eq. (|^). 
The zero-point fiuctuations are given by Eq. (^ with 
LUn replaced by cun- Of course, the harmonic approx- 
imation breaks down for the fundamental mode as we 
approach its critical strain. 



Ic - Iq _ /TTK^ 



(9) 



At critical strain, the effective potential for the funda- 
mental mode is purely quartic (see Fig. |l|) whereas the 
higher modes remain harmonic in leading order, with the 
first harmonic frequency being reduced by about 13% 
from its uncompressed value. 

To go beyond the harmonic approximation and calcu- 
late the excitation energies and quantum fluctuations of 
the fundamental mode, we consider the Hamiltonian 



H 



1 

2fl 



dx U^ + Viyix)], 



(10) 



where Il{x,t) = 6L/Sy(x,t) — iiy{x,t) is the canonical 
momentum. In the subcritical compression regime, a 
normal mode expansion of H leads to a phonon descrip- 
tion with interactions arising from the anharmonic term. 
These interactions occur physically because phonons 
stretch the bar. Even the zero-point motion has a stretch- 
ing effect, but this can be absorbed in the length pa- 
rameter I. Thus, at temperatures below the first har- 
monic threshhold, kT < ?i(I'2, the anharmonic effect on 
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the fundamental is its own self-interaction. So the effec- 
tive Hamiltonian for the fundamental mode, obtained by 
taking the ground state expectation value in all higher 
modes, is a quartic function of the "fundamental dis- 
placement" Y, the Fourier component of the fundamen- 
tal mode jl3|. The quantum states and associated en- 
ergy spectrum of fundamental vibrational states are then 
given by the Schrodinger equation, 



92 



2m* dY^ 
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^{Y)^E'^{Y) (11) 



where Em and '^^{Y) are the energy eigenvalues and 
eigenfunctions, and —ifid/dY = P is the momentum op- 
erator canonically conjugate to Y. 

The potential energy has the form of a Ginzburg- 
Landau free energy |]l^ , with strain playing the role of 
temperature: 



Er - e 



(12) 



The displacement Y is analogous to the order parame- 
ter, in that its classical equilibrium value vanishes below 
critical strain but takes a nonzero value Y — > ±Ymin — 
zb-y/jal//? above it, breaking the reflection symmetry of 
the Hamiltonian. Of course the quantum mechanical 
ground state has (Y) = 0, but sufficiently far into the 
double-well regime, this ground state is a superposition 
of macroscopically distinct states. Thus, monitoring the 
position of the bar on a time scale less than the tunnel- 
ing time would yield results clustered about just one of 
the potential minima, Ymm or —Ymin, not both. With 
this in mind, wc plot the ground and first excited state 
energies in Fig. ||, and in Fig. ^ the ground state quan- 
tum fluctuations Ay, as functions of the strain near its 
critical value. These are plotted in dimcnsionlcss energy 
and length units, E/Ec and AY/AYc, where Ec and AYc 
are the ground state values at critical strain, 



/ ^2 \ 2/3 / \ 

E, = 0A2{ — ) and AFe = 0.68 — - 

\m* J \m*p J 

and 7 is the departure from critical strain, 

1 ■ 



2 \ 1/6 



(13) 



£c \^Yj ' 



(14) 



scaled so that 7 = — 1 when the barrier height Vo = 
0^/4/3 is equal to Ei. In defining AY, care should be 
taken to account for the onset of a nonzero mean dis- 
placement as one moves into the double well regime. In 
the subcritical and near-critical regimes, the usual defi- 
nition AY^ = (Y"^) — (y)2 surely applies, where brackets 
indicate ground state expectation values and (Y) — 0. 
However, the comments above suggest that deeper in the 
double well regime, a more sensible practical definition is 



AY"^ = {min{Y ± YminY), the rms departure from the 
nearest potential minimum. Fig. ^ shows that quantum 
fluctuations defined in this way become much less than 
the well separation at fairly modest (negative) values of 
7. In the region 7 ~ — 1, the latter definition of AY"^ 
loses its physical meaning and larger fluctuations may 




FIG. 2: Ground and first excited state energies, Ei and E2, 
of the fundamental vibrational mode as functions of the de- 
parture from critical strain 7. The dotted curve shows the 
barrier height dependence on 7. 




FIG. 3: Curve A is the ground state fiuctuation Ay, calcu- 
lated using the full quartic potential. The solid portions of 
the curve are obtained from the calculations; the dot-dash re- 
gion is a guide to the eye. For comparison, curve B shows Ay 
obtained in the harmonic approximation. The dotted curve 
shows the position of the potential minimum, Ymin- 
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TABLE I: Table of excitation energies AE = E2 — E\ and rms midpoint fluctuations ISY for Si bars (linear dimensions I, d, 
w), and C nanotubes of length I and outer (inner) diameters d2 (di). The Young's modulus and density of Si are Q — 130 GPa 
and p = 5000 kg/m^. The values taken for C nanotubes are |^| Q = 1.8 TPa and p = 2150 kg/m"^. Energies AE are given in 
temperature units (note that IGHz — 48mK). 



Si bar Si bar C nanotube C nanotube 



;(nm) 500 50 /(nm) 500 50 

d(nm) 10 5 d2(nm) 10 5 

w(nm) 20 10 di(nm) 5 1 



relaxed 

A£;o(mK) 4.4 222 24.4 1100 

Ayo(A) .0061 .0054 .0093 .0082 



critical 

A£;c(mK) .027 2.0 0.21 14.5 

AYc{A) .055 .041 .073 .051 



be expected, as suggested by the dot-dash line on Fig. ^ 
which is only a guide to the eye. 

A useful perspective is provided by comparing the ex- 
act quantum fluctuations with those given by the har- 
monic approximation. Indeed, making this approxima- 
tion in the double well regime requires that one of the 
wells is ignored while fluctuations are defined about the 
minimum of the other, which clearly will become invalid 
as Ay approaches Ymm in the near-critical regime. As 
shown in Fig. ^, this approximation predicts that the 
fluctuation AY diverges as |e — ec\~^^'^- It is nonetheless 
a very good approximation outside a small region around 
the critical point (roughly —6 < 7 < 3), where the diver- 
gence is prevented by the quartic term in the potential 
energy. 

To address the magnitude of quantum fluctuations in 
real systems. Table | lists the first excitation energies 
AE = E2 — El and ground state (quantum) fluctua- 
tions Ay for rectangular silicon bars and cylindrical 
multiwalled carbon nanotubes. Numbers are given for 
two cases - the critical and the relaxed (uncompressed) 
states. The remaining Schrodinger equation parameters 
are calculated using hinged boundary conditions, with 
the results 



to/2 



and 



/3 = to*(wi/2k)2 



(15) 



The larger tubes have dimensions typical of the multi- 
walled tubes whose vibrational properties were studied 
by Treacy et al. The smaller tube is probably the 

smallest that would support buckling and retain its elas- 
tic integrity [T^ . Nanotube frequencies were found using 
Eq. (|) with = {dl + dl)/lQ, where di and ^2 are in- 
ner and outer diameters. In the critical and double well 
regimes, because of cylindrical symmetry, the nanotube 
states will be described by a Mexican hat potential rather 



with y replaced by a two-component vector. In order 
to make a comparison to silicon bars, the nanotube en- 
tries in Table | refer to a single Cartesian component of 
y. Fig. prefers specifically to one-dimensional transverse 
motion; the Mexican hat potential would have its main 
effect in the negative 7 regime where the first excitation 
energy E2 — Ei would fall less rapidly because the exci- 
tations are rotational rather than tunneling in character. 

For all systems shown, the zero-point fluctuations pre- 
dicted in the relaxed states are enhanced by nearly an 
order of magnitude by applying critical compression. On 
the other hand, the excitation energies are all reduced 
by about two orders of magnitude, so that thermal fluc- 
tuations can swamp the quantum contributions. On the 
other hand, the large differences between the relaxed and 
compressed energy scales suggests that the relaxed sys- 
tems could be supercooled by compression toward their 
critical points. For example, the smaller carbon nan- 
otube could be prepared initially very close to its ground 
state by cooling to 500 mK. Critical compression with- 
out heat transfer would then cool the tube to a few mK, 
at the same time enhancing its zero-point motion by a 
factor of about 6 (see Table ||) . Subsequent equilibration 
to 500 mK would then bring the tube to a "classical" 
equilibrium state with thermal fluctuation Ayt « 0.25 A 
|l7| , a further factor of 5 enhancement. 

Finally, we comment on the possibility of observing 
tunneling in nanomechanical systems such as those con- 
sidered here. We can roughly estimate the number of 
bound states with energy below the top of the barrier; 
it is of the order of the barrier height Vo divided by the 
level spacing TiCji. From Eqs. 0, n2 and we find 



N 



£c 



3/2 



(16) 



than a double well. In this case, Eq. I\n^ still applies. Taking the smaller bar, if we go to twice the critical com- 
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prcssion, e = 2ec, then iV ^ 3 x 10^. In order to tune 
the potential to hold about 10 bound states in each well, 
one would have to apply strain with extreme delicacy, 
e — Ec ^ 10~'*ec- Controlling the strain to this precision 
for sufficient time to identify tunneling, as distinct from 
thermal or other noise, will be difficult. Thus, while the 
observation of tunneling will likely be very challenging, 
the prospect of exploring tunable quantum fluctuations 
in this system, and the connection to Ginzburg-Landau 
theory, are intriguing. 
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